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We argue that (TMTSF) 2 PFs compound under pressure is 
likely a triplet superconductor with a vector order parameter 
d(k) ee (d a (k) / 0, d c (k) =?, d v Qt) = 0); |d a (k)| > |dc(k)|. 
It corresponds to an anisotropic spin susceptibility at T — 0: 
XV = Xo, Xa Xo, where xo is its value in a metallic phase. 
[The spin quantization axis, z, is parallel to a so-called b'- 
axis]. We show that the suggested order parameter explains 
why the upper critical field along the b'-axis exceeds all para- 
magnetic limiting fields, including that for a nonuniform su- 
perconducting state, whereas the upper critical field along the 
a-axis (a _L b') is limited by the Pauli paramagnetic effects 
[I. J. Lee, M. J. Naughton, G. M. Danner and P. M. Chaikin, 
Phys. Rev. Lett. 78,3555 (1997)]. The triplet order parame- 
ter is in agreement with the recent Knight shift measurements 
by I. J. Lee et al. as well as with the early results on a destruc- 
tion of superconductivity by nonmagnetic impurities and on 
the absence of the Hebel-Slichter peak in the NMR relaxation 
rate. 

PACS numbers: 74.70.Kn, 74.20.-z, 74.60.Ec 



Quasi-one-dimensional (Q1D) organic compounds 
(TMTSF) 2 X (X = PF 6 , C10 4 , etc.) have been intensively 
investigated since the discovery of superconductivity 1,2 in the 
first organic superconductor (TMTSF)2PF6. From the be- 
ginning, it was clear that their properties are unusual. It was 
found 3 " 8 that superconductivity in (TMTSF) 2 X (X = PF 6 , 
CIO4) is destroyed by nonmagnetic impurities. This was in- 
terpreted in terms of a possible triplet pairing of electrons 9 . 
Another unusual feature, the absence of the Hebel-Slichter 
peak in the 1/Ti NMR data in (TMTSF) 2 X (X = PF 6 , 
CIO4) 10, 11,12 , was prescribed 13 to the existence of zeros of 
a superconducting order parameter on the Q1D Fermi sur- 
faces (FS). As was stressed 13 , the early experiments 3-8 ' 10 ' 11 
provided information only about an orbital part of the order 
parameter and could not distinguish between some triplet and 
singlet pairings 2 ' 13 . 

To reveal triplet superconductivity, experimental tests 
which probe a spin part of an order parameter are essen- 
tial. Among them, are: a surviving of triplet supercon- 
ductivity in Q1D case 14-17 at magnetic fields higher than 
both the upper orbital critical field and the Clogston para- 
magnetic limit 18 , observation of spin- wave exitations 15 , the 
Knight shift measurements 12 and some others. Nowadays, 
interest in a possible triplet pairing has been renewed due to 
remarkable measurements of the upper critical fields (which 
are sensitive to a spin part of the order parameter) in 
(TMTSF) 2 C10 4 and in (TMTSF) 2 PF 6 at P ~ 6 kbar by 
Naughton, Lee, Chaikin and Danner 19 " 21 and due to the the- 
oretical analysis 16 of these experiments. The experimental 



fields along b'-axis (which are 3 times bigger 20,21 than the 
Clogston paramagnetic limit) were shown 16 to be even bigger 
than the paramagnetic limit 16 ' 22 for the Larkin-Ovchinnikov- 
Fulde-Ferrell (LOFF) phase 23 . Therefore, measurements 19-21 
were interpreted 16 ' 19-21 in term of triplet superconductivity. 
Recently, Lee et al. 12 have found no change of the Knight shift 
for H || b' in a superconducting phase of (TMTSF) 2 PF 6 at 
P ~ 6 kbar. This is consistent with the results 16 ' 19-21 and 
strongly supports the triplet scenario 9,16,19-21 of supercon- 
ductivity. 

The goals of our paper are: 1) To calculate the para- 
magnetic limited field along b'-axis, Hp , for the LOFF 
phase in a Q1D superconductor, taking account of both 
the paramagnetic 16 and orbital destructive effects against 
superconductivity. (We show that the calculated value of 
Hp is 4-5 times less than the experimental fields 20,21 in 
(TMTSF) 2 PF 6 ); 2) To demonstrate that the value of H b ' 
becomes consistent with 20,21 if we switch off the paramag- 
netic effects. (These indicate that an electron spin suscep- 
tibility along b'-axis, xv , at T = is equal to its value in 
a metallic state, xo> which is a distinct feature of triplet 
superconductivity 24,27 ); 3) To stress that the experimental 
critical fields 20,21 along the conducting chains (i.e., along a- 
axis), Hp, are strongly paramagnetically limited and thus the 
corresponding electron spin susceptibility Xa -C Xo at T — 0; 
4) To show that the above described properties are naturally 
explained within the framework of a triplet superconductiv- 
ity scenario with the following vector order parameter frozen 
into the crystalline lattice (i.e., the case of strong spin-orbit 
coupling 27 ): 

d(k) = (d a (k) jt 0,d c (k) =7,d b ,(k) = 0); |d (k)| > |d c (k)| 

(1) 

corresponding to the BCS-pair's wave function 
*(k) = [-da(k) + zd c (k)] I T T >+[da(k) + zd c (k)] I I I ) 

(2) 

and to the anisotropic spin susceptibility at T = 0: 

XV = Xo , Xa < Xo , (3) 

where [ "f \ (| 4 )) stands for a spin-up (spin-down) electron 
with respect to the quantization axis z || b' [a(x) _L b'(z) _L 
c*(y)], the momentum k defines the position on the FS. [We 
stress that b' is the easy axis for a spin direction in a spin- 
density-wave (SDW) phase of (TMTSF) 2 PF 6 . Thus, one may 
expect that the order parameter (1) is the most stable since it 
corresponds to the BCS pairs (2) only with Sy ss S z = ±1]. 
At the end of the paper, we discuss some consequences of a 
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group theory classification of the possible triplet phases, in- 
cluding the most probable orbital part of the order parameter 
and a possibility to break the time reversal symmetry. 

Q1D electron spectrum corresponds to two open sheets of 
the FS 1 ' 2 : 



e ± (p) = ±v F (p a TPf) - 2t b cos(p b b') - 2t c cos(p c c* 



(4) 



where +(— ) stands for the right (left) sheet of the FS; vf = 
t a a/y2 and pp are the Fermi velocity and Fermi momentum, 
respectively; t a ~ 1600 K, t b ~ 200 K and t c 5 K ; (h = 1). 

Singlet (S = 0) and triplet (S = 1) phases are characterized 
by the following wave functions of the BCS pair's 27 : 

^(k,r) = (|TI>-|IT>)^(k I r), S = ; (5) 

<Mk,r) = I TT )[-d»(k,r) + irf a (k,r)] + (| |I ) + | IT )) 

d z (k,r) + |U)[d x (k ) r)+id„(k,r)], S=l. (6) 

[In Eqs. (5,6), S 1 is the total spin of the BCS pair, r is 
its coordinate of a center of masses; tp s (k,r) — ip s (— k,r), 
d(k,r) = -d(-k 1 r)]. 

At H 0, ^(k, r) and d(k, r) do not depend on r. Electron 
spin susceptibility tensor, Xi,h at T = for a singlet phase is 
Xij = whereas for a triplet phase is given by 27 : 



X*,j = Xo ( 8i,j 



d?(k)d,-(k) \ 
d*(k)d(k) i 



(7) 



where dij = 1 if i = j and 8ij — if i 7^ j; <^|d(k)j 2 ^ k = 1, 
(•••), means an averaging over the FS. [Here, we consider only 
unitary triplet phases 27 (i.e., d a (k)d*(k) = d* (k)d c (k))]. 

At first we consider the case H || b'(z). In singlet phase 
(5), superconductivity is destroyed by paramagnetic effects in 
arbitrary directed magnetic field. In a triplet phase (6), as it 
follows from Eq. (7) , d b i (k) = d z (k) component is responsible 
for the deviation of the spin susceptibility xv = Xzz from xo- 
If dy (k) 7^ there exist two related phenomena: the paramag- 
netic destructive mechanism against superconductivity and a 
change of the Knight shift at T < T C (H). Let us calculate the 
upper critical field for H || b'. By using a common approach 28 
to the upper critical field of a clean superconductor 25 with 
open electron orbits and with one-component order parame- 
ter, it is possible to prove that Eq. (5) of Ref. [28] , 
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is extended to a singlet phase ip 3 {k., r) = /(k)A(x) as 
well as to the triplet phases di(k,r) = [d a = l,d c = 
0,d b , = 0)/(k)A(x) and d 2 (k,r) = (d a = 0,d c = 0,d b , = 
l)/(k)A(a;). [Here, <^|/(k)| 2 ^ k = 1; g is an effective electron 
interaction constant, d is a cutoff distance; a = ^/2t b /t a , 
ijj c = evFHc* /c, A = Atc/idc, \ib is a Borh magneton, e 
and c are the electron charge and the velocity of light, cor- 
respondingly; S z ~ 1 for singlet and for d2-triplet phases 



whereas S z — for di-triplet phase. By solving Eq. (8) 
numerically for S z = 1, a = 0.17, \dH b ' /dT\ Tc ^ 2 T/K, 
v F = 10 7 cm/sec, t c ~ 3K, T c (0) = 1.14 K, c* = 13.6 A (see 
Refs. [1,2,19-21,29]), we found that the calculated value of 
the paramagnetic limited critical field, Hp ~ 1.3 — 1.4 T, is 
4-5 times less than the experimental ones 20-21 (see Fig. ^]). 
A similar analysis for di-triplet phase (which is not param- 
agnetically limited) shows that superconductivity survives at 
H b ~ 6 T and T ~ 0.2 — 0.25 K in qualitative agreement with 
experiments 20 ' 21 (see Fig. |l|). On the basis of the calculation 
of Hi and H b , we can conclude that |d;/(k)| = |d z (k)| ~ 
in Eq. (7) and thus xv = Xzz — Xo- Note that the recent 
Knight shift measurements 12 are also in favor of xv = Xo 
below T C (H). 

If we consider the case H || a(x) then d a = d x - component 
of the order parameter (6) is responsible for the destructive 
paramagnetic effects against superconductivity and for the 
change of the Knight shift at T < T C (H) (see Eq. (7)). Let 
us calculate the critical field for H || a in di(r) = (d a 7^ 
0, d c 7^ 0,d b r = 0) A(x) triplet phase (which is paramagneti- 
cally limited). The corresponding linearized gap equation can 
be obtained from the common Eq. (5) of Ref. [28] : 
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where 7 = t a a/ (2t b b). Numerical solution of Eq. (9) (with 
the same values of parameters as Eq. (8)) shows that the 
best fitting of the data 20 ' 21 at H < 1.5 T (see Fig. 0) corre- 
sponds to S z — 0.9 (i.e., d a ^ 0.9, Xa — 0.2 xo <C Xo) and 
\dH a / dT\ Tc — 8 T/K. The latter is in a good agreement with 
the experimental slopes 20 ' 21 \dH b ' /dT\ Tc ^ 2 T/K since the 
value of tb/ta — 8.5 is known 29 . Note that the accuracy of 
our calculations does not allow us to distinguish between the 
triplet phases with d c — and \d a \ > \d c \. 

Summarizing, our analysis of the experimental critical 
fields 20 ' 21 measured in (TMTSF) 2 PF 6 at P ~ 6 kbar has 
shown that paramagnetic destructive effects against supercon- 
ductivity do not affect H b whereas H a is paramagnetically 
limited at H < 1.5 T. These are naturally explained within a 
triplet scenario of superconductivity 9 ' 16,19-21 with the triplet 
order parameter (1). We suggest to measure the Knight shift 
along the a-axis at H < 1.5 T and T < T C (H) to prove the or- 
der parameter (1). Note that temperature dependence of the 
critical field along a-axis, H a (T), changes drastically 20 ' 21 at 
H > 1.5 T. We speculate that at H > 1.5 T there may appear 
a triplet phase with d(k) _L H, which minimizes the magnetic 
contribution to the free energy 30 . Nevertheless, we cannot 
completely exclude another possibility - the appearance of the 
LOFF state at H > 1.5 T for H || a. Note that our theoretical 
analysis of the critical fields is based on the the Fermi-liquid 
picture 29 proved at P ~ 6 kbar in (TMTSF) 2 PF 6 . At higher 
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pressures, P ~ 9.8 kbar, the behavior of (TMTSF) 2 PF 6 may 
deviate from the Fermi liquid one 31 . 

At the end of the paper, we would like to make a few 
comments based on symmetry arguments. We classify the 
possible triplet phases in the case of strong spin-orbit cou- 
pling for orthorhombic (D2h) and triclinic (Ci) point group 
symmetries (see Table [i]) , where the matrix order parameter 
A(k) = di(k)fi, (fi = ioi<j y \ Oi are the Pauli matrices). As 
it seen from Table |, there are no degenerated orbital states, 
thus a time reversal symmetry is broken only if a nonunitary 
triplet phase appears 27 . In our particular case, this happens 
when d a (k)d*(k) 7^ d* a (k)a! c (k). Using the expression for a 
gap in a quasi-particle spectrum 27 , 8(k) = |d(k)| (the unitary 
case), it is possible to make sure that there are no generic 
phases with the lines of zeros on the FS in accordance with a 
common theorem 32 . This is in agreement with the experimen- 
tal data 26 ' 33 which seem to be in favor of fully gapped FS and 
against the existence of isolated zeros on theFS 32 . Therefore, 
we speculate that the orbital part of the order parameter is 
likely d a (k) ~ d c (k) ~ sgn(k a ) which corresponds to a fully 
gapped Q1D sheets of the FS. From Table |l[ it is possible to 
conclude that, for a triclinic space group of (TMTSF)2PF6, 
the most generic case is d a 7^ 0, d c 7^ and dy 7^ 0. How- 
ever, it is known 1,2 ' 34 that the spin dependent interactions in 
a SDW phase of (TMTSF) 2 PF 6 (which has a common bound- 
ary with the superconducting phase) result in an alignment of 
spins along b'-axis. Therefore, it is natural to expect the form 
(1) for the superconducting order parameter corresponding to 
the absence of the BCS pairs with Sy = (see Eq. (2)). 
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M. J. Naughton, K. Oshima, M. Sigrist, V. M. Yakovenko for 
useful discussions. A.G.L. is especially thankful to S. Brown, 
P. M. Chaikin, I. J. Lee and M. J. Naughton for fruitful and 
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TABLE I. Triplet order parameter A(k) for D^h, and d 
groups. (A~I are constants) 
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FIG. 1. Circles stand for the critical magnetic fields along 
b'-axis: open circles show experimental curve 20 ' 21 , a full circle 
corresponds to the calculated paramagnetically limited value 
of Hp at T = in a singlet superconductor whereas crossed 
circles show the calculated non-paramagnetically limited crit- 
ical fields H b (T) for a triplet order parameter (1). Trian- 
gles stand for the experimental critical fields 20 ' 21 along a-axis, 
Hp(T), (a _L b'). In the inset, the experimental values 20 ' 21 of 
Hp(T) are shown in comparison with the calculated paramag- 
netically limited field (full line) for a triplet order parameter 
(1) (see the text). 



